POLYNOMIAL INTERPOLATION ON THE UNIT SPHERE II 
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Abstract. The problem of interpolation at (n + l)^ points on the unit sphere 
S'^ by spherical polynomials of degree at most n is proved to have a unique 
solution for several sets of points. The points are located on a number of circles 
on the sphere with even number of points on each circle. The proof is based 
on a method of factorization of polynomials. 



1. Introduction 

— Jb \^ ti 2 



Let §^ = {a; : — 1} denote the unit sphere of M"^, where 



Let n„(S^) denote the space of spherical polynomials of degree n, which is the 
restriction of polynomials of degree n in three variables to It is known that 

dimn„(§2) = (n + 1)^ n > 0. 

The problem of interpolation on the unit sphere by polynomials is as follows: 

Problem 1. Let X = {sLi : 1 < i < (n + 1)^} be a set of pairwise distinct points 
on S'^. Find conditions on X such that there is a unique polynomial T £ n„(§^) 
satisfying 

T(a,)-/„ a,eX, l<i<(n+l)2, 
where {fi} is an arbitrary set of data. 

If there is a unique solution to the interpolation problem, we say that the problem 
is poised and that X solves Problem 1. This problem has been studied recently in 

OS El mini HI Hi 

Although almost all choices of X will solve Problem 1, it is difficult to know 
whether a given set X will work since computing the determinant of the interpola- 
tion matrix is difficult. In (14) a large family of sets of interpolation points is given 
explicitly, each set solving Problem 1. Let us briefly describe this construction. 
The (n + 1)^ points lie on n + 1 distinct latitudes (parallel circles on §^), and each 
latitude contains an odd number of equidistant points. The number of points needs 
not to be the same on each latitude and there is no restriction on the position of 
the latitudes. For the simplest case n = 2m, the set of {2m + 1)^ points lie on 
2m + 1 latitudes, each of them containing 2m + 1 equally spaced nodes. In 
another family of points that solves Problem 1 was found, for which n = 2m — 1. 
There the points lie on 2m latitudes and each latitude has an even number of 2to 
equally spaced points. In this case, the 2to latitudes are divided into two groups; 
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the equidistant points on one half of the latitudes need to differ by a rotation from 
the points on the other half of the latitudes. While the proof in ^ is based on 
the analysis of the determinants of the interpolation matrix, the proof in |14) uses 
a factorization method which avoids the determinants. Furthermore, the factor- 
ization method provides many more sets of points leading to poised problems. A 
key observation in 14 is that the use of equidistant points allows us to reduce the 
problem on the sphere to a special trigonometric interpolation problem. 

The purpose of this paper is to show that the factorization method also works 
in the setting of an even number of points on each latitude. Again, the use of 
equidistant points reduces the problem to an interpolation problem of one variable. 
However, the new interpolation problem is different from the one with an odd 
number of points on each latitude and has to be solved using a completely different 
method. In comparison to 4 , the factorization method allows to obtain more sets 
of points that solve Problem 1. 

For the background of polynomial interpolation in general, we refer to the survey 
article and the references therein, even though interpolation on the sphere is not 
discussed there. Let us also mention that the factorization method is closely related 
to the method used for polynomial interpolation on the unit disk in ^EllHl- Apart 
from a result in Jj,, which is a simple consequence of Bezout's theorem, the family 
of points found in 0] , and those stated below appear to be the only ones that 
are given explicitly for all n. 

The paper is organized as follows. The factorization method is studied in Section 
2 and its application to polynomial interpolation on is given in Section 3. 

2. Factorization of polynomials 

2.1. Polynomial representation. For fixed a G (—1, 1), let §^(a) :— {{x,y,z) G 
: z — a} denote the circle on resulting from the intersection of with the 

plane z — a. This set is called latitude at z = a. 

On the unit sphere §^ it is more convenient to work with spherical coordinates, 

X — sin 9 sin (j>, y ~ sin 6 cos 0, z — cos 6, < (p < 2tt, < 9 < tt. 

For a polynomial r„ G n„(§^), we introduce the notation T„ defined by 

T„(6',(?!)) = T„(sin6l COS0, sin6l sin0, cos6l), < </)< 27r, O<0<tt. 

li X = {{xi,yi, Zi) : 1 < i < M} is a set of points on we also use the notation 
X = {{6i, cf)i) : 1 < i < M} for the corresponding set of spherical coordinates. 
It has been shown in Section 2 of (14) that the polynomial T„ can be written as 

(2.1) f„(0,0) =ao(cos^) 

n 

+ [aA;(cos6') (sin6')''cosfc(/) + 6fc(cos6')(sin6')''sinfc(/)] , 
fc=i 

where afe(-) and bk{-) are polynomials of degree n—k in one variable. Note that for 
any fixed 0, the polynomial T„(9, •) is a trigonometric polynomial of degree n. 

Below we will consider interpolation problems based on points that are equidis- 
tantly distributed on an even number of latitudes, each of them containing an 
even number of nodes. To describe these points, it is convenient to introduce the 
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following notation: 

e„,. := : = M±^, j=0,l,...,2.s-l 

where s G N and a G [0,2). These points can be considered as equidistant points 
on the unit circle using the mapping (j) i— > e*"^. The parameter a indicates that the 
points are defined up to a rotation by an angle of a-n /2s. 

Lemma 2.1. Let n = 2m— 1 and a G [0, 2). For (j) G Qa,m, 

(2.2) f„(0,0) =ao(cos0) 



m — 1 

('r.^o fl^/'c.^,. , fr,r.c a\fc;^ a\'^"i-k\ 

k=l 

+ (6fc (cos 6*) (sin 9)'' + V2m-k{cose) (sin Q)^"'-'' ) sin kcj)] 



+ 22 [(afc(cos6')(sin6')''+W2m_fc(cos6')(sin6')^™"'') cosfci^ 

1 



+ \ am{cos9) cos bm{cos9) sin — j (sin 6')™ cos ym4> ^ y , 

where, for fc = l,...,m— 1, 

U2m-k{t) = a2m-kit) COS ttTT + b2m~k{t) siu aTT, 

V2m-k{t) = a2m_fe(t)sinQ;7r- 52m-fc(f) cos a7r 
are polynomials of degree fc — 1 . 

Proof. We split the sum in H2.1|l into two sums, one over 1 < fc < m — 1 and the 
other over m < k < 2m — 1. In the second sum we change the summation index 
k t-^ 2m— k and use the elementary relations 

cos(2m — k)(t> = cos((2j + 0)11 — kef) — cos an cos fc</) + sin an sin k(j>, 

sin(2TO — k)(f> — sin((2j + a)n — k(f>) — sin an cos k(j) — cos an sin k4>, 

which holds for (f> G Qa,m- Combining the two sums, we obtain 



T„(0,0) = ao(cos6l) + V [(afc(cos6l)(sin6i)'' + U2m-fe(cos6i)(sin6i)2™-'=)cosfc(/) 



m— 1 

r/„_ fl\/„;„ ^3^fe I „. fl\/„:„ ^3^2m-fc^ 

k=l 

+ {bk (cos 61) (sin 9)''+V2rn-k (cos 61) (sin 6»)2"-'= ) sin k(f>] 



+ {a,n{cos9) cos(Q!7r — m(f>) + 6„i(cos0) sin(a7r — m(j))){sin9y 

Using the addition formula for the cosine and the sine function 

(an\ an . / an\ . an 
mcj) ) cos h sin I m(j) ) sin — , 

m0 1 cos cos I m0 ) sin — , 

2 / 2 V 2 / 2 

the Ora and &,„ terms of the above expression of T!„ can be rewritten as 



^am(cos0) cos 6m (cos 0) sin (sin 6*)™ cos ^rrif/) — 

(an an \ / 

a,„ (cos 0) sin — +6„i(cos6') cos— j (sin 0)™ sin (^m0 - 



) 

an\ 



Consequently, formula 1)2. 2|l follows from the fact that (p G &a,m satisfies sin(m0 — 
an/2) = 0. ' □ 
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Lemma 2.2. LetT2,n-i be given as in (|2.2(l and 9 e (0,7r). If T2m-i{0 , (j)) = for 
(j) S 0Q,m, then ao{cos6) = 0, and 

ak{cosd) + (sin6')^''"^^'^(a2m-fc(cos6') cosaTr + 62m-fc(cos0) sinaTr) = 0, 

bk{cos9) + (sin6')^'"^^''(a2m-fe(cos6') sinaTr — b2m-k{cos6) cosaTr) = 0, 

for 1 < fc < m — 1 . Furthermore, 

ceTT air 

(2.4) am(cos6') cos Om(cos6') sin = 0. 

Proof. The proof uses the following fact. Interpolation on the 2m zeros of the 
function sin {jn(f> — ^) inside [0, 27r) by a trigonometric polynomial of the form 

rn— 1 

Sm-i{4>) = io + (flfc COS fc6' + bk sin fc6') + a„ cos(m(/) — a7r/2) 
fe=i 

is unique (see ^1 Vol. II]). Observe that the points in Qa,m are exactly zeros of 
sin (rn0 — an /2) and cos(m0 — an /2) = ±1 for e O^.m. 

By Lemma im T2m-i takes the form of 12.2|l . The assumption that T2m-i{6, </>)= 
implies that the coefficients of r2m-i(^, •) are all zero. This leads to ao(cos6') — 0, 

afc(cos6')(sin6l)'= +U2„-fc(cos6l)(sin6l)2™^'= = 0, 
&fc(cos0)(sin0)^ + V2m-fe(cos0)(sin0)2'"-^ = 0, 

for k — l,...,r7i — 1, which become the equations in (|2.3|l upon multiplying by 
(sin 6')^'^, and 

/ an anx 

(^Om (cos 0) cos 6m(cos 0) sin j (sin0)™ =0, 

which gives ((2.4|l . Note that sin0 as £ (0,7r). □ 

To proceed from here, we want to choose 2to distinct such that whenever the 
polynomials in H2.3|l and H2.4|l vanish on these 2m points, they will be identically 
zero. To this end, however, we need to impose an additional symmetry. We choose 
di to satisfy 

(2.5) 6'2,„+i-i TT - e*,, G (0, tt), 1 < i < m. 

In other words, we choose the latitudes to be symmetric with respect to the equator. 

If p(t) is a polynomial of degree n, we denote by and 'p°'^^ the even and the 
odd part of p, respectively. To be precise, if p(t) ~ X]j=o Qj^''? then 

p™™(t)= a2,e' and p°'^^(0 = ^ a^-it^'-^ . 

0<2j<N l<2j-l<N 

Lemma 2.3. Let T2m-i be given as in 1)2.2(1 . If for some 9 in (0, tt) 

T2m-i{0, 0) = 0, e eo,m and T2,n-l{n -0,(j>) ^0, cj) <^ 6l,m, 
then, setting t = cos9, we have ao{t) — and, for 1 < fc < m — 1, 

w/iere either p2„i-k-i{t) = ak{t) and qk^i{t) = a2m-k{t), or p2m-k-i{t) = bk(t) 



an- 



d qk-i{t) = -b2m-k{t)- Furthermore, a.m{t) = bm(t) = 0. 
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(2.7) 



Proof. The assumption allows us to use the previous lemma. Since a = or a = 1, 
the fact that am{t) — bm{t) = follows immediately from l|2.4(l . For 1 < fc < to — 1 
and a = 0, the equations (|2.3|l become 

ak{cos6) + (sin0)^™"^''a2„_fc(cos6') = 0, 

fefc(cos(?) - (sin^^)2™-2fc62,„_,.(cos^?) = 0. 

For l<fc<TO— = l and 6 replaced by n — 6, the equations (|2.3|l take the 
form 

afe(-cos(?) - (sin^^)2™-2'=a2™-fc(-cos0) = 0, 



(2.8) 



Since p{t) +p(— t) = 2p''™"(t) and p{t) — p{—t) = 2p°'*^(t), combining equations 
(j2.7|l and p. 8(1 proves the result. □ 

We use the notation P2m-k-i and qk-i since they are polynomials of degree 
2to — fc — 1 and k—1, respectively. In the following we will work with the equations 
in 1(2. 6|l . If p is an even polynomial, it can be written as p{t) = p*{t^); if p is an odd 
polynomial, it can be written as p{t) = tp*{t^). Thus, in place of ((2. (ill we will need 
to consider polynomials of the form p{t) + tq{t){l — t^Y and tp{t) + q{t){l — t^Y . 
We will need to study the possibility of interpolation by such polynomials. This is 
discussed in the following subsection. 

2.2. Chebyshev systems. A family of functions . . . , (t>r} is called a Cheby- 
shev system on a set i? C M, if every linear combination from the span{(/)i, . . . ,(j)Y\ 
has at most r zeros in E] in other words, interpolation on r points by functions in 
the span{(/)i, . . . ^(pr} has a unique solution. In this subsection we prove that the 
families of functions in ((2.6(1 are Chebyshev systems on (0, 1). 

Proposition 2.4. Let r and s he two nonnegative integers such that r > s > 0. 
For e — or 1, let 

g{t)=pr{t')+t^\l-t'Y-'qs-i+e{t^), 

where p^ and qg^m^^ are polynomials of degree r and ,s — 1 + e, respectively. If g 
vanishes on r + s + 1 + e distinct points in (0, 1), then g{t) = 0. 

Proof. We first prove the case that the power of in g{t) is taken as t and e = 0. 
The cases t~^ or e = 1 are similar; in fact, the proof for the case e = 1 is identical, 
and only minor changes (merely the numbers below will change) are needed for 
the case that t^^ is taken as t~^. 

Changing variables t^t^ shows that we need to prove that if 

g*{t)^Pr{t) + Vt{l-tY-'qs-l{t) 
vanishes on r + s + 1 distinct points in [0, 1], then g*{t) = 0. Let 

Using Rolle's theorem repeatedly, we see that it suffices to prove that if h{t) vanishes 
on s distinct points in (0, 1), then qs-i{t) = 0. 
Since qs-i is a polynomial, we can write it as 

qs-i{t) =bo + bi{l-t) + ... + 6,_i(l - tY'^. 
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Using the Leibnitz rule repeatedly, we have 



r — k—j 



In the following, we will use the convention that Jli=a ~ ^ whenever b < a. This 
leads to 



h{t) 



s~l r—s 



fc=0 3=0 



3 



k=0 



where 



h* 

Ui, 



k — 1 s — k—1 r — s 

— [|(2z + l) Y[ (2i-l) and hk{t) ^Y.^'^'J*' 



+ k 



i=0 i=l 

in which the coefficients a^j are given by 



(2.9) 



r — s 

i 



\[{2k + 22 + 1) J|(2(r -k-i)-l). 



i=0 



'■=3 



We note that all coefficients akj are positive numbers. The polynomial hk is of 
degree r — s + fc. In order to prove the proposition, we need to show that the set 
{/iQ, hi, . . . , hg^i} forms a Chebyshev system on (0, 1). In other words, we need 
to prove that the matrix (^j(ifc))j7-=o invertible for any set of distinct points in 
(0,1). 

Let t — {ti,t2, ■ ■ ■ ,ts} be a given set of distinct numbers in (0, 1). For a given 
set of nonnegative integers A = {jo,jii ■ ■ ■ ijs-i}, we introduce the notation 



V{X;t) = dct 



4-30 
'■2 



■^1 ''2 



In the case of A = {s — 1, s — 2, . . . , 0}, we denote the determinant by Vs{t), which 
is the Vandermonde determinant 

T4(t)=det[ti+J^-.Lo= n fe-*^)- 

l<i<j<s 

For a given set of nonnegative integers A, we further introduce the notation 

V{X;t) 



SA(t) 



Vs{t) 



Note that s\ is a symmetric polynomial in t and s\ is zero if jo,ji, ■ ■ ■ , Js-i are 
not pairwise distinct. If /i = (/^oiMi • • • iMs-i) is a partition, that is, Mo ^ Mi ^ 
■ ■ ■ ^ Ms-i 0, Mi G N, and ji — fi^ + n — i + I for < j < s — 1, then sx is 
called a Schur polynomial, cf. |21 for details. It is known that Schur polynomials 
can be written as a linear combination of monomial symmetric polynomials and 
the coefficients in the linear combination are all positive (called Kosta numbers). 
For our purpose, it is enough to note that the Schur polynomials are positive when 
t; > for all 1 < Z < s. In particular, it follows that if jo < ji < . . . < js-i, then 
s\{t) is positive when ti > for all 1 < ^ < s. 
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Using the definition of the polynomials hj we can write 
(2.10) ^^detih,{tk+i)r.;,l^ 

r — s r — s+1 r — 1 

= X! "' X! aOjo°^lji-l • ■ •as-lJ.-l-(s-l)SJ0Jl,■■■J=-l(*)• 

J0=0 Jl=l Ja_i=S-l 

For A = {jo, ji, ■ • • , js-i} we will also denote the coefficient of s\ in the above sum 
by A\; that is 

— aojo«i,ii-i ■ • ■ Q^s-ija_i-(s-i)- 
Recall that Uj^k are all positive numbers. If A = {jo, ji, • ■ ■ j js-i} is not a partition, 
then a proper permutation of jo, ji, • • ■ , js-i will be. The determinant changes sign 
when two rows are exchanged, so that s\ is positive if the permutation is even and 
it is negative if the permutation is odd. Every permutation can be factored into a 
number of transpositions. A transposition means exchanging two elements. If sx 
is negative, there is a transposition of A, call it A', such that sx'{t) = —sx{t) > 0. 

Let Sx (t) be negative and assume that jp and jg are a pair in A = { jo , ji , • ■ • , js- 1 } 
such that p < q but jp > jg. Considering the summation indices in (|2.1U|I . we must 
have jp > jq > q > p. Let A' be the image of A under the transposition (p, q) , that 
is, with jp and jg exchanged. Then sx' (t) > 0. The coefficients of these two terms 
in (|2.1U|) . Ax and Ax' differ by two terms only. We have 

Ax' ~ Ax — {o.p.jq-pO'q.j-p-q ~ ^p,jp-p^q,jq-q) W O.i.ji-i- 

i¥'P;q 

We now show that Ay > Ax, which will complete the proof of the proposition. 
Recall the definition of ak,j in H2.9|l . Let us denote bkj — (ikj / {^^'') ■ Then it is 
easy to verify that 

bp,j,-P ^ (2g-l)(2g-3)...(2p+l) ^ 

bq,j^-q (2s-2g-l)(2s-2g-3)...(2s-2p+l) ' 

which is independent of jp and jg. Consequently, we have 



^''^\jp-q+^ J'"\ jp-P 'J a^jp-? 
which implies that Ax' > Ax and completes the proof. □ 

Proposition 2.5. Let m and k be integers such that 1 < k < m. Let p2m-k-i cind 
qk-i be polynomials of degree 2m — k — 1 and degree /c — 1 , respectively. Lfti, . . . ,tm 
be distinct numbers in (0, 1) and 

(2.11) PT^i'Lk^iiu) + qfMm t^T-^ = 0, 

(2.12) P2^-^-k{^^) + €-l{U)(\ - t^r-' = 0, 
then p2m^k-i{t) = and qk-i{t) = 0. 

Proof. Depending on k being even or odd, we need to consider the following four 
cases. 
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Case 1. A: is even. Setting r — m — {k + 2)/2 and s = {k — 2)/2, equation ^2.11^ 
becomes 

Prit^) + t, qs{tl){l - tjy-' =0, 1 < i < r + ,s + 2. 
From ProDOsition l2.4l with e = 1 and t^^ — t, it follows that Pr{t) = and qs{t) = 0. 

Case 2. A: is even. Settine; r = m — {k + 2)/2 and s = {k - 2)/2, equation 
becomes 

Uplitl) + g:(t2)(l - tjy-^ = 0, 1 < I < r + 5 + 2. 

Multiplying the equation by we can use Proposition 12 .41 with e = 1 and t^^ = 
to conclude that Pr{t) = and qs(fy = 0. 

Case 3. k is odd. Setting r = m - (/c + l)/2 and s = (fc - l)/2, equation (|2.11(l 
becomes 

vAtl) + i,g:_i(t,2)(l - i^-' - 0, 1 < z < r + s + 1. 
From Proposition l2.4l with e = and = t, it follows that Pr{t) = and qs{t) = 0. 

Case 4. A: is odd. Setting r = m - {k + 2)/2 and s = (A: - 2)/2, equation (|2.12() 
becomes 

+ q*s-iitf){l - t?)''"^ = 0, 1 < z < r + 5 + 1. 

Multiplying the equation by t^^ , we can use Proposition l2.4l with e = 1 and t^^ — t 
to conclude that Pr{t) = and qs{t) = 0. □ 

2.3. Factorization method. The following factorization theorem holds the key 
to our main result. 

Theorem 2.6. Let m and s be positive integers satisfying m < s < 2m— 1. Denote 
X — s—m+1. Let 01, . . . , 02A be distinct numbers in (0,7r) such that 92x+i~i — Tr—Oi 
for j = 1, . . . , A. Denote 

X = {{9i, 4>ij) : 4>i^j e Qo.mi 1 < i < A, and e Qi.m, X + l <i < 2A}. 

Lf Ts G ns(§'^) satisfies 

then there is a spherical polynomial T'*_2A ^ ns_2A(§^) such that 

T{x,y,z) = - cosei),T*^2\{x,y,z). 

1=1 

In particular, T*_2x = if s — 2m—l. 

Proof. We start with the formula (|2.1(l . which becomes 

s 

fn{0,(j)) =ao(cos6')+^ [afc(cos6') (sin6')'' cosfc0 + 6fe(cos6')(sin0)'' sinfc^] , 
fc=i 

where afe(-) and 6fe(-) are polynomials of degree s — k. For i = 1, 2 . . . , A, we can 
follow the proof of Lemma [2. 21 and Lemma f2. 31 and distinguish the following three 
cases. 

Case 1. For < A; < 2™ - s - 1, 

afc(cos6',)(sin6',)'' = 0, 5fc(cos 61,) (sin 6',)'' = 0, i = l,...,2A, 
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setting 60 = 0. Since 9i e (0, tt), we have that ak{cos9i) — and bk{cos9i) = 0, for 
i = 1, . . . , 2A. Reeah that ak and bk are polynomials of degree s — k > s — (2m — 
s — 1) = 2A — 1. Consequently, there exist polynomials and 6^, both of degree 
s — k — 2A, such that 

2A 2A 

ak{t) = Y[{t- cos 9,)al{t) and bk{t) ^Y[{t - cos 9,)bl{t). 

i=l 1=1 

In the extreme case k = 2m — s—1, we have a2m-s-i — b2m-s-i — 0. 
Case 2. For k = m, 

{sin9i)"^am{cos9i) — Q and i — 1,...,A, 
(sin6'i)'"6„(cos6l,) = and i A + 1, . . . , 2A. 

Since 9i G (0, tt) and both a,„ and 6„ are polynomials of degree s — m = A — 1 
vanishing at at least A points, they have to be identically zero. 

Case 3. For 2m — s<fc<m— 1, we end up with equations similar to 1)2. 6|l : 

p°'^i{u) + gir-(.-fe)-2(iO(i - 1^-^-^^+' = 0, 

for i = 1, 2, . . . , A, where either Ps-k{t) = ak{t) and q2X-{s~k)-2{t) = a2m~k{t), or 
Ps^k{t) = &fe(0 and q2A-(s-fc)-2(^) = -&2m-fc(<)- In deriving the above equations 
we have used several times the identity A = s — m + 1. Recall that and b^ are 
polynomials of degree s — k; the subscript of the polynomials p2A-fc and q2\-(s-k)-2 
again indicate their degree. 

It is easy to see that the system of equations (|2.13ll is exactly the one being stud- 
ied in the previous subsection, namely H2.11|l and H2.12|l . Hence, using Proposition 
12 .51 we conclude that 

ak{t) = and a2m-k(t) =0, k ~ 2m — s, . . . , m — 1, 
bk{t) = Q and &2m-fc(i) = 0, fc 2m - s, . . . , m - 1. 

Together, these three cases show that we have the factorization 

2A 

fs{9,^)=\{{cos9-cos9{) 

1=1 

(2m-s-2 \ 
ag(cos6') + ^ (a^ (cos 0) COS fc(?i + 6^ (cos 0) sin fc0) j , 

which completes the proof. □ 

Using factorization repeatedly, we can obtain a complete factorization of a poly- 
nomial of degree 2m — 1 in n2m_i(S^). 

Theorem 2.7. Let n be an odd positive integer, a €z N, and Ai, . . . , Ao- be positive 
integers. Define Uk ~ Uk-i — 2Afc, for I < k < a, with ng = n. Assume that > 
forl<k<a-l. IfTne n„(§2) satisfies 

T„(0.,fc, (i)^,J,k) =0, 1 < z < 2Afc, < j < 2(nfc_i - Afc + 1) - 1, 1 < fc < a, 
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where 9i^k are pairwise distinct angles in (0,7r) with 9i,2\^+i-i = tt — 9i,i, I = 
l,...,Afe, <f)ij^k e 6o,nfc_i-Afc+i for 1 < i < Xk and £ Bi,„^_i-Afc+i for 

'^fc + 1 < * < 2Ai:, then there exists a polynomial T*^ G n„^(§^) such that 

Tn{x,y,z) = Yi Yii^ - cos6',,fc)T*^(a;,?/,2:). 

k=l i=l 

In particular, Tn{x, y, z) = if < 0. 

Proof. We apply the factorization result in Theorem 12 . 61 repeatedly with s ~ nk~i, 
m = rifc-i — Afc + 1 and A = Xk for fc = 1, 2, . . . , cr. □ 

Just as in the case of an odd number of points on each latitude (cf. JJ), the 
interpolation nodes in the above theorem are located on latitudes split up into a 
groups {S^(zi,fc) : 1 < « < 2Afc}, I < k < a, Zi^k = cos6'i,fc, and Z2Afc+i-i,fc = -Zi,k- 
Latitudes in different groups contain a different number of nodes. More precisely, 
each of the latitudes in the fc-th group, S^(zi^/c), S^(z2,fc), ■ • ■ ,^'^{z2\k,k), contains an 
even number of 2(nk-i — Xk + 1) equidistant points and the points lie on symmetric 
latitudes. In other words, points on S^(z2Afc+i-i,fc) and §^(zi,fc), « = 1, . . . , Afe, differ 
by a rotation of an angle of 7r/(2(nfc_i — Afc + 1)). 



3. Interpolation on the sphere 

Our main result on interpolation follows from the factorization Theorem l2.7l The 
following formula can be used to verify that the number of interpolation conditions 
matches the dimension of the polynomial space: 

(3.1) n^(S^) =dimn^_2A(§^) + 2A(2s-2A + 2). 

Theorem 3.1. Let n be an odd natural number and let Xi^ . . . , Xa (c G N) be 
positive integers, such that 

(3.2) Ai + ... + A, = 

Define Uk — Uk-i — 2Xk, for 1 < fc < tr — 1, with uq = n. Let 

X = {(6l,,fc, fc) : 1 < i < 2Afc, < J < 2(nfc_i - Afc + 1) - 1, 1 < fc < cr}, 

where 9i_k, 1 < J < 2Afc and 1 < k < a, are distinct numbers in (0,7r) with 
92\k+i-iM = TT - 9i^k (i = 1, . . . , Afc), 0jj^fc e 6o,nfc_i-Afc+i, for 1 < i < Xk, and 
4'i.j,k G ©i.Mfc-i-Afc+i, for Xk + l < i < 2Afc. Then the set X solves the interpolation 
problem in n„(§^). 

Proof. First, we verify that the dimension of n„(§^) matches the number of inter- 
polation conditions. Let \X\ denote the number of points in AC. It follows from 
equation H3.1|l that 

\X\ =^2Afc (2nfc_i-2Afc + 2) 

k=l 
a 

= ^(dimn„,_,(S2) - dimn„,_,_2A.(§')) = [n + lf = dimn„(S2). 

k=l 
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Thus, it is sufficient to show that if Tn G n„(S^) vanishes on X, then Tn{x, y, z) = 0. 
Under the condition %\.2\ . it follows that 

cr-l 

= n^-i - 2\„ = n„^2 - 2Acr-i - 2Ao- = • • • = n - 2 ^ Afc = -1 < 0. 

fe=i 

Hence, the factorization in Theorem 12 . 71 shows that T{x, y, z) = 0. □ 

For a fixed n — 2m— 1, this theorem contains a number of different interpolation 
processes. In fact, for each positive odd integer n, the number of point sets X 
which can be deduced from Theorem 13 . II depends on the partition number of (n + 
l)/2. Every solution of equation H3.2|l leads to a set of points defining a poised 
interpolation problem in n„(§^). The number of solutions of such an equation 
grows exponentially as n goes to infinity. Moreover, the order of Ai, . . . , matters; 
i.e. , different permutations of a solution Ai, . . . , A^- of equation H3.2(l give different 
sets of interpolation points. 

Among the solutions of equation 1)3. 2|l . one extreme case is it = 1, for which the 
equation has only one solution Ai = (n + 1)/2. In this case, the interpolation points 
are located on ri + 1 symmetric latitudes §^(2:1), §^(^2), ■ • ■ ,S^(-Zn+i), each of them 
containing n + 1 equidistant points. This case has already appeared in 0] Theorem 
2.5]. 

The other extreme case is a — {n + l)/2 and Ai = • • • = ^{n+i)/2 = 1- There, 
Uk — n ~ 2k and we have {n + l)/2 groups of two symmetric latitudes, where the 
ones in the kth group contain 2(n — 2fc + 2) equidistant nodes. The points on a 
latitude are rotated by an angle of 7r/(2(n — 2k + 2)) with respect to the points on 
the corresponding symmetric latitude of the same group. 

Example 3.2. To illustrate the power of the factorization method, we present the 
possible point distributions for n ~ 3,5 and 7. 

• 71 = 3 

(1) cr = 1, Ai = 2; 4 latitudes each with 4 points, 

(2) cr = 2, Ai = A2 = 1.' 2 latitudes with 6 points and 2 latitudes with 2 
points. 

• n — 5 

(1) cr = 1, Ai = 3; 6 latitudes, each with 6 points. 

(2) a ~ 2, Xi + X2 = 3 has two solutions. 

(a) Ai = 1,A2 — 2; 4 latitudes with 8 points and 2 latitudes with 2 
points. 

(b) Ai = 2, A2 = 1.' 2 latitudes with 10 points and 4 latitudes with 4 
points. 

(3) cr = 3, Ai = A2 = A3 = 1; 2 latitudes with 10 points, 2 latitudes with 6 
points and 2 latitudes with 2 points. 

• n = 1 

(1) cr = 1, Ai = 4; 8 latitudes with 8 points. 

(2) cr — 2, Ai + A2 = 4 has three solutions, 

(a) Ai = 2, A2 = 2; 4 latitudes with 12 points and 4 latitudes with 4 
points; 

(b) Ai = 1, A2 = 3; 2 latitudes with 14 points and 6 latitudes with 6 
points; 
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(c) Ai = 3, A2 = 1.' 6 latitudes with 10 points and 1 latitudes with 2 
points. 

(3) cr = 3, Ai + A2 + A3 = 4 has three solutions, 

(a) Ai = 1,A2 = IjAs = 2; 2 latitudes with 14 points, 2 latitudes 
with 10 points and 4 latitudes with 4 points. 

(b) Ai = 1,A2 = 2, A3 = 1; 2 latitudes with 14 points, 4 latitudes 
with 8 points and 2 latitudes with 2 points; 

(c) Ai = 2,A2 = 1,A3 = 1; 4 latitudes with 12 points, 2 latitudes 
with 6 points and 2 latitudes with 2 points. 

(4) cr = 4, Ai = A2 = A3 = A4 = 1; 2 latitudes with 14 points, 2 latitudes 
with 10 points, 2 latitudes with 6 points and 2 latitudes with 4 points. 

It is well-known that interpolating polynomials can be used to construct cuba- 
ture formulas on the unit sphere (cf. jl2|). In fact, integrating the interpolation 
polynomial in n„(§^) yields a cubature formula on the sphere which is exact for 
spherical polynomials of degree ri. Among the point sets in Theorem 13. II the case 
where the points are distributed on 2m symmetric latitudes, with each latitude con- 
taining 2m equidistant points, is of particular interest. In this case, the cubature 
formula is simple and can be explicitly given. 

Proposition 3.3. Let m be a positive integer. Let 61, ... , 62m be pairwise distinct 
numbers in (0, tt) with 02m+i-i = it — 9i, i ^ 1, . . . ,m, and a £ {0, 1}. Then for 

aiiT2m-i e n2„-i(§2), 

„ m2m — 1 2m 2m — 1 

/ T2,„_i(0dw(0 = -E^^ E ^2m-l(0.>°) + - V A,; V Tsm-l (0. , ) 
•'^ i=l j=Q i=m+l j=Q 

where (j)" = (2j + a)Tr/2m, and Ai is given by 

/I 2m „ 
-r-r t — COS Ok 

-Ifcii^. COS0.-COS0fc 

Proof. Let the interpolation polynomial r2m-i be of the form (|2.1|l . We use the 
quadrature formula 

1 p2-ir 1 2m-l 

which is known to hold for every trigonometric polynomial of degree at most m 
(see, for example, [EI Vol.2, p. 8]). Using formula 1)3. 3|l and the interpolation 
property of T2m-i, it follows that 

1 ^27r _ ^ 2m-l ^ 

ao(cOS0,) - — / T2m-l{0^,(|))d^ ^ ^ T2m-M,<t>'^), 

j=0 

for every fixed 9i, 1 < i < 2m. Consequently, oq, which is a polynomial of degree 
2m — 1 in one variable, is uniquely determined by these 2m interpolation conditions. 
It follows that 



m j 2m- 1 \ 2m / 2m- 1 \ 

«0W - E I ^ E ^2m-l(e..<^°) I m+ r2m-l(^.0l) | ^t). 



J=0 
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where £i{t) = Jlfc^i fc/i(^ ~ cos 6k)/ {cos 9i — cos 6k)- Using the change of variables 

T2„,-mduj{0= / f2,n-i{6,(p)sin6d(pd6, 
Jo Jo 

the integral of T'2m-i over the surface of the sphere is equal to 

/ T2™_i($)dw(0 = 27r / ao(cos6) sin 9d9 ^2tt ao{t)dt. 
JS2 Jo J~i 

The stated formula follows from the formula for ao{t) given above. □ 

In particular, this result shows that the cubature formula is nonnegative, if cos 6i 
are chosen so that are nonnegative. This holds, for example, if cos^i are the 
zeros of the Legendre polynomial P2m of degree 2m, or the zeros of a quasi Legendre 
orthogonal polynomial P2m + af2m-i with mild conditions imposed on a G M (see, 
for example, ^S])- In QDji the positivity of the cubature in this case has been 
proved by working directly with the interpolation matrix. 
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